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Indian Standard 

METHODS FOR 
FIXING THE CAPACITIES OF RESERVOIRS 

PART IV FLOOD STORAGE 
0. FOREWORD 

0.1 This Indian Standard ( Part IV ) was adopted by the Indian Standards 
Institution on 19 February 1971, after the draft finalized by the Catchment 
Area and Reservoirs Sectional Committee had been approved by the CivH 
Engineering Division Council, 

0.2 Flood storage depends on the height at which the maximum water 
level (MWL) is fixed above the normal conservation level (NCL). The 
determination of the MWL involves the routing of the design flood through 
the reservoir and spillway. When the spillway capacity provided is low, the 
flood storage required for moderating a particular flood will be large and 
vice versa, A higher MWI involves larger submergence and hence this 
aspect has also to be kept in view while fixing the MWL and the flood 
storage capacity of the reservoir. 

©♦3 This standard consists of four parts and the other parts are as follows: 

IS: 5477 (Part I )-1969 Methods for fixing the capacities of reservoirs: 
Part I General requirements 

IS:5477(Part II)-1969 Methods for fixing the capacities of reser- 
voirs: Part II Dead storage 

IS: 5477 (Part III)-1969 Methods for fixing the capacities of reser- 
voirs ; Part III Live storage 

0*4 For the purpose of deciding whether a particular requirement of this 
standard is complied with, the final value, observed or calculated, express- 
ing the result of a test or analy^ shall be rounded' off in accordance with 
IS : 2-1960*. The number of si :tuficant places retained in the rounded off 
value should be the same as thai of the specified value in this standard. 



1. SCOPE 

1.1 This standard ( Part IV ) covers the criteria and procedure to be follow- 
ed in fixing the flood storage capacity of a reservoir consistent with the 



•Rules for rounding off numerical values ( revised ). 
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safety of the structure itself and the life and properties downstteam of the 
reservoir. 

2. TERMINOLOGY 

2.0 For the purpose of this standard, the following definitions shah 
apply. 

2.1 Normal Conservation Level (N<3L) — The normal conservation 
level is the highest level of the reservoir at v^rhich water is intended to be 
held for various uses other than flood control, 

2.2 Full Reservoir Level (FRL) — Is the highest level of the reservoir at 
which water is intended to be held for various uses including part or total 
of the flood storage without allowing any passage of water through the 
spillway. 

2.3 Maximum Water Level ( MWL ) — Is the highest level to which the 
reservoir waters will rise while passing the design flood with the spillway 
facilities in full operation, 

2.4 Surchavge Storage — The storage between the crest of an uncontroll- 
ed spillway or the top of the crest gates in normal closed position and the 
MWL is termed as the surcharge storage. 

2.5 Maximum Probable Flood — Is the flood that may be expected from 
the most severe combination of critical meteorologic and hydrologic 
conditions that are reasonably possible in the region, and is computed by 
using the maximum probable storm which is an estimate of the physical 
upper limit to storm rainfall over the basin. This is obtained from storm 
studies of all the storms that have occurred over the region and maximizing 
them for the most critical atmospheric conditions {see A-2AA of Appen- 
dix A also). 

2.6 Standard Project Flood — Is the flood that may be expected from the 
most severe combination of meceorologic and hydrologic conditions consi- 
dered reasonably characteristic of the region and is computed from the 
standard project storm rainfall reasonably capable of occurrence over the 
basin in question and may be taken as the largest storm which has occurred 
in the region of the basin during the period of weather record. It is not 
maximized for most critical atmospheric conditions but it may be transposed 
from an adjacent region to the water-shed under consideration. 

2.7 Design Flood — Is the flood adopted for design purposes. It may be 
the maximum probable flood or the standard project flood or a flood 
corresponding to some desired frequency of occurrence depending upon the 
standard of security that should be provided against possible failure of the 
structure. 
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3. GENERAL 

3.1 A prerequisite essential to the safety of the dam is a decision in regard 
to the standard of security that should be provided against possible failure 
of the concerned structure during extraordinary floods. In case, failure of 
the structure is likely to result in loss of life or widespread property damage 
in addition to the damage of the structure itself as in the case of spillways 
of large reservoirs, a very high degree of security shall be provided against 
failure under the most severe flood conditions considered reasonably 
possible. But when loss of life is not involved, economic considerations 
would prevail to govern the selection of the design flood. 

3.2 It should be assumed that the reservoir would be filled to the full 
reservoir level at the beginning of the spillway design flood. Even though 
the contemplated plan of reservoir operation indicates that a portion of the 
storage capacity below the FRL probably would be available at the 
beginning of the spillway design flood, the possibility of improper operation 
of regulating outlets as the result of incorrect flood predictions, mechanical 
difficulties, plugging of conduits by debris, or negligent attendance 
may justify the assumption of a full reservoir at the beginning of the design 
flood. Moreover, future developments may require revisions in the original 
plan of reservoir operation or changes in the use of the reservoir that would 
increase the probability of a full reservoir at the beginning of the spillway 
design flood. 

4. ESTIMATION OF DESIGN FLOOD 

4.1 The methods in vogue for estimation of design flood are broadly 
classified as under: 

a) Application of a suitable factor of safety to maximum observed 
flood or maximum historical flood, 

b) Empirical flood formulae, 

c) Envelope curves, 

d) Frequency analysis, and 

e) Rational method of derivation of design flood from storm studies 
and application of unit hydrograph principle. 

4*2 Application of a Suitable Factor of Safety to Maximum Observed 
Flood or Maximum Historical Flood — The design flood is obtained by 
applying a safety factor which depends upon the judgement of the designer 
to the observed or estimated maximum historical flood at the project site or 
nearby site on the same stream. This method is lirnited by the highly 
subjective selection of a safety factor and the length of available stream flow 
record which may give a quite inadequate sample of flood magnitudes likely 
to occur over a long period of time. 
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4.3 Fmpirical Flood Formulae' — The empirical formulae commonly 
used in the country are the Dicken's formula, Ryve's formula and Inglis' 
formula in which the peak flow is given as a function of the catchment area 
and a coefficient. The values of the coefEcient vary within rather wide. 
limits and have to be selected on the basis of judgement. They have 
limited regional application, should be used with caution and only, when a 
more accurate method cannot be applied for lack of data. 

4»4 Envelope Carves — In the envelope curve method maximum flood is 
obtained from the envelope curve of all the observed maximum floods for a 
number of catchments in a homogeneous meteorological region plotted 
against drainage area. This method, although useful for generalizing the 
limits of floods actually experienced in the region under consideration, 
cannot be relied upon for estimating maximum probable floods for the 
determination of spillway capacity except as an aid to ludgement. 

4,5 Frequency Analysis — ^The frequency method involves the statistical 
analysis of observed data of a fairly long (at least 25 years ) period. 

4«5*1 A purely statistical approach when applied to derive design floods 
for long recurrence intervals several times larger than the data has many 
limitations and hence this method has to be used with caution. 

4*6 Rational Method of Derivation of Design Flood from Storni 
Studies and Application t>f Unit Hydrograph Principle 

4.6.1 The steps involved, in brief, are: 

a) analysis of rainfall Vs run-off* data for derivation of loss rates under 
critical conditions; 

b) derivation of unit hydrograph by analysis ( or by synthesis, in cases 
where data are not available ) ; 

c) derivation of the design storm; and 

d) derivation of design flood from the design storm by the application 
of the rainfall excess increments to the unit hydrograph. 

Appendix A provides the criteria for estimation of design flood. 

4.6.2 Unit Hydrograph — Limitations 

a) The unit hydrograph principle is not applicable for drainage basins 
having an area of more than 5000 km^ where valley storage effects 
are not reflected and where variation of rainfall in space and time 
shows a tendency to become too great to be reflected in the unit 
hydrograph. 

b) Application of the unit hydrograph principle is also not recom- 
mended for catchments having an area less than about 25 km^. 
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c) Large number of raingauges suitably located should be available in, 
the entire catchment to reflect the true weighted rainfall of the 
catchment. 

d) Unit hydrograph principle is not applicable when appreciable 
proportions of the precipitation occurs in the form of snow or wh^ 
snow covers a significant part of the catchment. 

5. DETERMINATION OF MAXIMUM WAIER LEVEL 

5»1 The maximum water level of a reservoir is obtained by routing the 
design flood through the reservoir and spillway. This process of computing 
the reservoir stages, storage volumes and outflow rates corresponding to a 
particular hydrograph of inflow is commonly referred to as flood routing. 
The routing is determined by the following: 

a) Initial reservoir stage ; 

b) The design-flood hydrograph ; 

c) Rate of outflow including the flow over the crest, through sluices or 
outlets and through power units ; and 

d) Incremental storage capacity (each stage). 

5.1.1 Although there are definite relations between reservoir inflow, 
storage, and outflow these relations are usually diflicult to express algebrai- 
cally. Therefore, a step-by-step computation procedure is followed 
whereby the increase in storage and rate of outflow resulting from the 
volumes of inflow during successive short increments of time are computed. 
Increments of inflow are computed for periods of time sufficiently short to 
warrant the assumption that the mean of the inflow rates and the mean of 
the outflow rates at the beginning and end of the intervals would closely 
approximate the average rates for the respective periods. 

5.2 Appendix B provides an illustration with the details of routing the in- 
flow hydrograph through a storage reservoir and fixing the maximum water 
level. Appendix G gives an example based on graphical ( Sorensen ) 
method of routing the inflow hydrograph through a reservoir. 



APPENDIX A 

( Clauses 2.5 and 4.6.1 ) 

CRITERIA FOR ESTIMATION OF DESIGN FLOODS FOR 
DAMS AND OTHER HYDRAULIC STRUCTURES 

A-1. GENERAL 

A-1.1 Each site is individual in its local conditions, and evaluation of 
causes and effects. Therefore, only general guidelines are provided and 
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the hydrologists and the designer would have the discretion to vary the 
criteria in special cases, where the same are justifiable on account of assess- 
able and acceptable local conditions, which should be recorded and 
have the acceptance of the competent authority. 

A-2. DETERMINATION OF DESIGN FLOOD 

A-2.1 Major and Mediam Reservoirs 

Ao2.1«l Maximum Probable Flood — In the design of spillways for major and 
medium reservoirs (with storages more than 6009 hectare metres) the 
maximum probable flood should be used. The maximum probable flood 
is estimated from the maximum probable storm applying the unit hydro- 
graph principle. The maximum probable storm is an estimate of the 
physical upper limit to storm rainfall over a basin. It is obtained from 
the studies of all the storms that have occurred in the region and maximiz- 
ing them for possible moisture charge and for storm efficiency. 

The method of moisture adjustment commonly used involves the 
estimation of air moisture content from surface level dew point observa- 
tions. Maximizing storm efficiency is achieved by storm transposition 
which assumes that the most effective combination of ^torm efficiency and 
inflow wind has either occurred or has been closely approached in the 
outstanding storm on record. Further adjustment of storm precipitation 
data to the estimated maximum sustained wind to carry this maximum 
moisture supply into the project basin is not attempted unless very high 
design safety factors are desired or only very limited storm rainfall data 
are available. 

A-2«l«2 Distribution of Storm Rainfall — The distribution of storm intensi- 
ties for small durations is obtained on the basis of recorded data of self- 
recording raingauge stations in the catchment or region. 

A-2.1 .3 Maximization for Unit Hydrograph Peak — It is observed in practice 
that the unit hydrograph peak obtained from heavier rainfalls is about 25 to 
50 percent higher than those obtained from the smaller rainfalls. There- 
fore, the unit hydrograph from the observed floods may have to be suitably 
maximized up to a limit of 50 percent depending upon the judgement of the 
hydrologist. In case the unit hydrograph is derived from very large floods, 
then the increase may be of a very small order; if it is derived from low 
floods the increase may have to be substantial. 

A-2.1.4 Loss Rate — The loss rate should be estimated from the volume 
of observed flood run-off and the corresponding storm that caused the 
flood. A minimum loss rate should be worked out and adopted for com- 
puting the design flood. 

A-*2.2 The probability method, when applied to derive design floods for 
long recurrence intervals several times larger than the length of data, has 
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many limitations. In certain cases, however, like that of very large catch- 
ments where unit hydrograph method is not applicable and where suffi- 
cient long term discharge data is available, the frequency method may be 
the only course possible. In such cases the design flood to be adopted for 
major structures should have a frequency of not less than once in 1 000 
years. Where annual flood values of adequate length are available, they 
are to be analyzed by Gumbel's method {see Appendix D), and where the 
^ata is short, either partial duration method or regional frequency techni- 
qv^e is to be adopted as a tentative approach and the results verified and 
che'^ked by hydrological approaches. 

Sometimes when the flood data is inadequate, frequency analysis of 
recorded storms is made and the storm of a particular frequency applied 
to the unit hydrograph to derive the flood; this flood usually has a return 
period greater than that of the storm. 

A-2.3 Barrages and Minor Dams — In the case of parmanent barrages 
and minor dams with less than 6000 hectare metres storage, the standard 
project flood or a K)0 year flood, whichever is higher, is to be adopted. 

A-2.3.1 For pick-up weirs a flood 50-100 years frequency should be 
adopted according to its importance, and level conditions. 

A-3. PROVIS'ONS FOR OTHER FACTORS 

A-3.1 The initial reservoir level before the impact of the spillway design 
flood has to be taken as at full reservoir level. In regions experiencing 
prolonged floods where storms in quick succession are experienced, the 
spillway may also be checked for design flood preceded or succeeded by a 
flood of once in 25 years frequency. The interval between these two 
floods ( peak to peak ) may be taken as 3 or 5 days according to as the 
region lies in an annual rainfall zone of more than 100 cm or less than 
100 cm respectively. 

A-3,2 To provide for mechanical and other failures, it is assumed that some 
gates as inoperative with a maximum of 10 percent and a minimum of one 
gate. For this purpose the designer may be permitted to incjrease permis- 
sible stresses treating it as an extraordinary occurrence, like earthquake. 

APPENDIX B 

{Clause 5.2) 

STEP-BY-STEP METHOD OF RESERVOIR FLOOD ROUTING 
FOR FIXING MAXIMUM WATER LEVEL 
(MODIFIED PULS' METHOD) 
B-1. BASIC DATA REQUIRED 
B-1.1 The following is the basic data: 

a) Reservoir level versus spillway capacity (outflow) {see Fig. 1), 
and storage capacity {see Fig. 1 ); 
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b) Full reservoir level; 

c) Inflow hydrograph {see Fig. 2); 

d) Initial outflow; and 

e) Initial storage (or initial reservoir level). 
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Fig. 1 Reservoir Levels Vs Spillway and Storage Capacities 
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Time Vs Reservoir Levels for the Design Flood Hydrograph 
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B*2. BASIC EQUATION FOR ROUTING FLOOD THROUGH 
RESERVOIR 

B-2.1 Basic equation for routing flood through reservoir is given below: 

where 

/ = time interval, 
/i = initial inflow in m^/s, 
/g = final inflow in m^/s after t h, 
Oi = initial outflow in m^s, 
O2 = final outflow in m^/s after / h, 
5*2 = final gross storage in cumec h after t h, 
Si = initial gross storage in cumec h, and 
/\5= incremental storage in time t h. 

Storage capacity is usually given in hectaremetres which should be 
converted into cumec h or cumec day, if the time interval / is expressed in 
hours or a part of a day. 

B-3. PROCEDURE 

B-3.1 Taking known values on one side 

(i._±A)^(^_|)„(| + o,) <„ 

1 cumec day = 24 cumec h — 8*64 ha,m ... -• (3) 

5ha,m =2-771 5cumech ... ... ...(4) 

5 =2-771 5/^ cumec h ... ... ...(5) 

— ha.m ^ ' 



t 



Taking the time interval / as 6 h 
X =0-463 5 ... ... ... ... (6) 

B-3.2 From Fig. 1, Table 1 is prepared. 
B-3.3 By means of Table 1, Fig. 3 is prepared. 
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TABLE 1 


RESERVOIR LEVEL Vs 
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( ROUTING PERIOD - 6 HOURS ) 








( Clause B-3,2) 
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B-3.4 Referring to Table 2, col 1 and 2 are obtained from the given 
inflow hydrograph ( Fig. 2 ) . The values given in col 3 are the average 
of successive inflows. Column 4 is obtained by subtracting the initial 0/2 
value from the initial 5/^ value ( which is either known or read out from 
Fig. 1 knowing the initial reservoir level). 

Column 5 is the sum of col 3 and 4* 

Column 6 is obtained by referring to Fig, 3 for the values of col 5. 

Corresponding to the values of col 6, the values in col 7 are obtained 
by referring to Fig. 1 . 

With the values given in col 6 and 7, the outflow hydrograph and the 
curve of reservoir level with time are drawn in Fig. 2. 
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TABLE 2 ROUTING OF DESIGN FLOOD THROUGH RESERVOIR 
(ROUTING PERIOD. 6 HOURS) 



( Clause B-3.4 ) 
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APPENDIX C 

( Clause 5.2 ) 

GRAPHICAL (SORENSEN'S) METHOD OF RESERVOIR 

FLOOD ROUTING FOR FIXING MAXIMUM 

WATER LEVEL 

C-1. BASIC DATA REQUIRED 
C-Ll The following is the basic data: 

a) Reservoir level versus spillway capacity (outflow) {see Fig. 4), 
and storage capacity {see Fig, 4); 

b) Full reservoir level; 

c) Inflow hydrograph ( see Fig, 5A ); 

d) Initial outflow; and 

e) Initial storage (or initial reservoir level). 

C-2. BASIC EQUATION FOR ROUTING FLOOD THROUGH 
RESERVOIR 

C-2.1 Basic equation for routing flood through reservoir is given below: 

dS:={I—0) dt ... ... ... ...(7) 

where 

dS= change in reservoir storage, 
/— rate of inflow, 

= spillway discharge rate (or outflow rate\ and 
dt = time interval. 

As the inflow is a function of time and outflow is a function of storage, 
equation (7) may be written as: 

S = Ss-S,:^I-^^ + ^^^^C.i ...(8) 

where 7 is the average inflow for the period l\t and the subscripts 
H' and '/' refer to the initial and final conditions of the period l^t. 

C-3. PROCEDURE 

C-3.1 For graphical solution, equation (^8) is transposed as follows: 

(5.+ ^) + (7-a)A* = 5,+ 0,^l ...(9) 
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Fig. 4 Reservoir Level Vs Storage and Oijtflow Curves 



A constant value of A ^ is selected, depending on the shape of the 

inflow hydrograph ( here A < ~ 1 h ) . The values of 5 + ~|^ are plotted 

as a function of (a) resen/oir level as shown in Fig. 5B, and (b) spillway 
discharge (or outflow) as shown in Fig. 5C. 

NoTS — The scales for the ordinates of Tig, 5A and 5C should be the same. 

In Fig. 5C by the side of S + ^ versus outflow curve, a line A is 

C 
drawn with the same origin such that its slope is ■ . where C is the 

ratio of *S + "^ scale to the scale ( that is, ratio of the scale of X-axis 
to the Y-axis of Fig. 5G ). 
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For this particular example, ^ f is taken as one hour and the scale 
1 250 
ratio works out to -ggo"'^ ^' which determines the slope of the line A, 

At 0800 h the reservoir inflow and the outflow are known. To deter- 
mine the reservoir level one hour later, that is at 0900 h, the following 
procedure is adopted: 

a) Draw a horizontal line a from the value of the average inflow 7 
for the period 0800 h to 0900 h, towards Fig. 5C. 

b) Draw a line b in Fig. 5G parallel to line Ay from the intersection 
of horizontal ( towards Fig. 5G ) corresponding to the initial out- 
flow at 0800 h and the curve of 5 H -— versus in Fig. 5G. 

c) The intersection point of lines a and b is projected vertically in 
Fig. 5B and 5G to get the points Tand Z respectively. 

d) The point T of Fig. 5B is projected horizontally to Fig. 5D to get 
the reservoir elevation at 0900 h. 

e) The point X of Fig. 5C is projected horizontally to Fig, 5A to eet 
the outflow Of at 0900 h. 

f ) With new reservoir level and outflow values as initial conditions 
the ^teps a to e are repeated for the next hourly increment. Thus 
reservoir stage hydrograph Fig. 5D and the outflow hydrograph 
Fig. 5 A (superimposed on the inflow hydrograph) are obtained. 

The curves of Fig. 4 are not directly involved in the graphical solu- 
tion but are the basic data for plotting Fig. 5B and 5G {see Table 3). 



TABLE 3 


COMPUTATION OF VALUES FOR Fig. 
( WITH At AS I HOUR ) 


5B AND 5G 


Elhvation 


Storage, 
1 000 ha. 


m 


OtTTPLOW, 

lOOOha. rr 


i/h 


2 
1000 ha. 


in 


J 000 ha. m 


(I) 


(2) 




(3) 




w 




(5} 


240-70 


1-500 




0-032 




0016 




1-516 


242-25 


2-100 




0-120 




0-060 




2-160 


243*02 


2-500 




0-200 




0100 




2-600 


243-80 


3*237 




0-302 




0-151 




3-388 


245-35 


5- 188 




0-612 




0-306 




5-494 


246-90 


7-500 




0-987 




0-494 




7-994 
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5C Spillway Discharge Vs S + 
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Fig. 5 Graphical (Sorensen's) MExkoD for Fixing Maximum Water Level 
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APPENDIX D 

[Clause A'2.2) 

ESTIMATION OF PARAMETERS OF GUMBEL'S EXTREME 

VALUE T3ISTRIBUTION AND DETERMINATION OF PEAK 

FLOW OF LONG TERM RETURN PERIOD 

D-1. MAGNITUDE OF PEAK FLOW OF A GIVEN RETURN 
PERIOD 

D-1.1 Gumbel postulated the use of extreme value distribution 

F{x)dx =^ ae-^'^^-^h'^ dx ... ...(10) 

where a and u are to be evaluated from n years annual observed peak 

flow values x^^ x^^ Xz^ , ;c„. 

Peak flovy magnitude xt of 7" years return period is evaluated from 

xt = u — {\ja) \oge loge T/( T — 1 ) 

=. „ _ 2;322_6r Q.3g2 3 + logic logic 'f^A ^-(1 U 

D-2. EXAMPLE 

D-24 A worked out example for data of Sabarmati at Dharoi ( 1935-52 ), 
as given in Table 4 illustrates the computational procedure 
(I mVs= 1000 1/s). 

Referring to Table 4, 

mean flood =^= -^-^= IfZ?, == 98-352 9 ... .-(12) 

n 17 ^ ^ 

Standard deviation s is given by 

Ex^ - {E xY 

256002 '--164446:1116_ ^3^22-242 650 



" 16 

.-. J =:= 75-645 506 

Computation for * a ' 

As a first approximation to a, take 

ai^ ^ X —= 0-016955 ... ,..(14) 

Then col 5 in Table 4 is completed. For ^2, «3, ^4, etc, of col 6, 7 and 8 
of Table 4, a step-by-step method of computing successive approximations 
Ck^i from ak is given in Table 5, which is self- explanatory ^ Successive 
approximations are continued in Tables 4 and 5 till a negligibly small 
value approaching of Ajt is reached in Table 5. 

1^ 
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TABLE 4 VALUE OF 


,-«* 












{ Clause D-2A) 








Si. 


Yeab 


Peak 


j^ 


r**i* 


r**2* 


rv 


rv 


No. 




Flow 
















Ratb (*) 
(10000 

1/s) 












(I) 


(2) 


(3) 


(4) 


f5) 


(6) 


(7) 


(8) 


1 


1935 


76 


5 776 


0-275 7 


0-263 1 


0-2624 


0-2624 


2 


1936 


3 


9 


0-950 4 


0*948 6 


0-948 6 


0-948 6 


3 


1937 


125 


15 625 


0-120 I 


0-1113 


0M10 8 


0M10 8 


4 


1938 


79 


6 241 


0-262 


0-249 7 


0-249 


0-249 


5 


1939 


227 


51529 


0021 3 


0-018 6 


0-018 4 


00184 


6 


1940 


91 


8 281 


213 7 


0-202 2 


0-201 6 


0-201 5 


7 


1941 


42 


I 764 


0-490 6 


0-478 2 


0-477 4 


0*477 4 


8 


1942 


125 


15 625 


0-120 1 


O-Ul 3 


0-1108 


0-110 8 


9 


1943 


193 


37 249 


0-037 9 


0033 7 


0-033 5 


0033 5 


10 


1944 


246 


60 516 


0*015 4 


0-013 3 


0013 2 


0-013 2 


11 


1945 


187 


34J969 


0042 


0-037 5 


0*037 2 


0037 2 


12 


1946 


45 


2 025 


0-466 2 


453 6 


0-452 9 


0-452 9 


13 


1947 


57 


3 249 


0-380 5 


0-367 5 


0*366 7 


0*366 7 


14 


1948 


20 


400 


0-712 4 


. 0-703 7 


0-703 3 


0-703 3 


15 


1949 


48 


2 304 


0-443 2 


0-430 3 


0*429 6 


C-429 6 


16 


1951 


102 


10 404 


0-177 4 


0-166 6 


0-166 I 


0*166 I 


17 


1952 


6 


36 


0-903 3 


0-899 9 


0-899 7 


0*899 7 




TOTAt 


1672 


256002 


5-632 2 


5*489 1 


5 4812 


5-481 1 



Computation for *tt' 

From the steady value so reached of a^ (or k ^ 4 ( in the present 
case ) , compute 

1/- Ic^elOX -— - r 1(^10 («)-logio (r*"'**''^)! ,..(15) 

= 64-297 104 

Putting these values of u and the steady a^ in (11), expected peak values 
XT for r= 50; 100; 200; 500, etc, years return periods are evaluated as 

x&Q = 285-97; jckw = 325"63; ^^200 = 365-15; x^oa = 417-30 

In order to safeguard against incidental errors of estimation arising 

from small samples of 29, 30, ..' years observed data available the 

above best estimates may be increased by 1*645 SE{xt) to warrant 95 
percent dependability of the design estimates. The standard error of 

fl 
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estimation SE{xt) is given by 

SE{x_t) = ^Pl + -V( 1-0-577 216 

-log, ioger/(r-i)rj ... ...(16) 

= — ^ fl + 0-607 927 1 1 0-422 7M 

-log. log. r/(r- 1)1']* ... ...(17) 

Thus 

X60 - 2 859 700 l/s ( = 2 859*7 m'/s) + 797-3 m'/s 

^100 = 3 256 300 „ ( = 3 256*3 mVs) + 916-1 m'/s 

a:20o--3 651 500 „ ( = 3 651-5 mVs) + 1 035*7 m^/s 

ATsoo =s 4 173 DOO „ ( = 4 173-0 m»/s ) + 1 194-6 m^/s 

Values of ^= — ioge log. Tl{ T— 1 ) for the most commonly used 
values of T like 20, 50, 100, 200, 500,.and 1 000 years are reproduced below 
for convenience of users: 

r-20 50 100 200 500 1000 

^ = 2-9702 3-9020 4-6002 5*2958 6-2137 6-9073 

D>2.2 It is advisable to procure always at least 20-25 years observed records 
and not to compute estimates for return periods longer than about 20w or 
25« years. 

l>-2.3 The steps remain very much the same irrespective of the value of n, 
the number of years. Though the computational steps by successive 
approximation method appear very arduous, their schematically outlined 
procedure as in D-2.1, even with recorded data of around 50 years should 
not take a person equipped with an ordinary calculating machine and a 
book of logarithm tables more than 6-7 hours. With electronic fast 
computers the stabilized expression can be secured in two to three minutes 
for which programme is given in D-3. 

D-2.4 The calculation of a 500 years ' flood, for example does nat tell when 
the flood is coming; it might occur in any year within that period or not 
until 500 years have elapsed. 

D-2.5 The choice of a suitable return period on which to base the * design 
flood' is the engineer's responsibility, and rests ultimately on his judgement 
and experience. The size and cost of the dam, the design freeboard, the 
amount of water stored, and the likely consequences of failure are factors 
which influence the selection. 

21 



TABLE 5 FOK VALUES OF 0;^ ^ i FROM m„ 

( Clause D'Zl) 



Step 


Quantity 


Jt^l 


Jt = 2 


*-3 


*-* 


1 


flk 


016 955 


0017 565 


0-017 602 


0017 602 58 


2 


a* logioe 


0007 363 


0-007 628 


0*007 644 


0-007 645 


3 


l/«* 


58-979 6 


56-931 4 


56*811 7 


56-809 9 


4 


!/«•* 


3 478-593 2 


3 241-184 3 


3 227-569 3 


3 227-364 7 


5 


X - 1M» 


39-373 3 


41-421 5 


41-541 2 


41-543 


6 


2/' 


5-632 2 


5-489 1 


5-4812 


5-481 1 


7 




240-075 2 

■■ 


228-368 4 


227-7111 


227-702 


8 




19 898-861 8 


18 563-741 2 


18 488-027 9 


18 487*199 8 


9 


f M - (7) ^ (5) (6) 


18-316 900 


1*001644 


0-015 475 


0000 663 


10 


/'(«*) -(8) + (5) (7) 
- (4) (6) 


-30 038-441 5 


-26 895-564 3 


-26 719-588 4 


-26 717*284 3 


U 




0000 610 


0000 037 


0000 000 58 


0000 000 02 


12 


flk + 1 = ^ifc + ** 


017 565 


0-017 602 


0-017 602 58 


0-176 026 



! 



<0 
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D-3, COMPUTER PROGRAMME FOR PEAK FLOOD ESTIMA- 
TION BY MAXIMUM LIKELfflOOD METHOD 

FIRST DATA CARD IS NO. OF OBSERVATION 4COLS, NAMES OF RIVER, 

SITE, UNIT OF OBSERVATION 8COLS EACH FOLLOWED BY DATA CARDS 

4COLS EACH OBSERVATION. 

DIMENSION X( 100 ), EX( 100 ) 

SUM = 0. 

SQ.SUM=0. 

READ 2, N, RVR, SITE, UNIT 
2 FORMAT (14, 3A8) 

READ3, (X(I),I = 1,N) 
3 FORMAT (20F4.2) 

DO 5 1=1, N 

SUM==SUM+X{I) 
5 SQ,SUM=SQ,SUM+X(I)**2 

B=N 

SQ,SUM^(SCiSUM— SUM*SUM/B)/(B-1.) 

SUM=SUM/B 

T=3.141592654 

A=T/SQRTF(6.*SQ,SUM) 

H=0. 

DO 1.1=1,20 

SIG-0. 

XSIG=0. 

X2SIG=0. 

A-A— H 

DO 8 1=1, N 

EX(I)c=EXPF(-A*X(I)) 

SIG - SIG + EX( I ) 

XSIG = XSIG + X( I )*EX( I ) 
8 X2SIG = X2SIG + X( I )*X( I )*EX( I ) 

FA == XSIG - ( SUM - l./A )*SIG 

FH = ( SUM - l./A ;*XSIG - X2SIG - SIG/A**2 

H ^ FA/FH 

IF(ABSF(H)-l.E-06)30, 30, I 
1 CONTINUE 
30 PIE = 6./( T**2 ) 

T= l./A 

U = T*LOGF( B/SIG ) 

X20 = U + T*2.970I86 

X50 = U + T*3.901953 

XIOO^U + T*4.600150 

X200 = U + T*5.295775 

X500 = U + T*6.213675 

P- l./(A*SQ,RTF (B) ) 

SEX20 = P*SQRTF ( 1. + PIE*( 0.422784 + 2.970186 )**2 ) 

SEX50 = P*SQ.RTF ( 1. + PIE*( 0.422784 + 3.901953 )**2 ) 

SEXIOO - P*SQ,RTF ( 1. + PIE*( 0.422784 + 4.600150 )**2 ) 

SEX200 = P*SQ.RTF ( 1. + PIE*( 0.422784 + 5.295775 )**2 ) 

SEX5G0 = P*SQRTF ( 1. + PIE*{ 0.422784 + 6.213675 )**2 ) 

PRINT 310, UNIT, RVR, SITE 

{Continued) 
23 
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310 FO(<MAT(IH1/////30X 43HEXPECTED VALUES WITH VARIOUS RETURN PERlOl 
1 46X, IH*, A8, 4X, 1H*//40X, A8,3X,2HAT,A8 ) 

PRINT 309 
309 F0RMAT( //29X. 8H20 YEARS, 12X, 8H50 YEARS, IIX, 9HI00 YEARS, MX, 9H2( 
I EARS, IIX, 9H500 YEARS) 
PRINT 311, X20, X50, XIOO, X200, X500 

311 FORMAT{8X.9HESTIMATE-,5 (4X. Fl^.9)) 
PRINT 312, SEX20, SEX50, SEXIOO, SEX200. SEX500 

312 FORMAT( /9X, 8HSTD ERR ^, 5 ( 4X, F16.9 ) ) 
STOP 

END 
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INTElJEMATlONAL SYSTEM OF UWITS (SI UNITS J 



TJQia 

iE-mi;KfrAture 
LyiQinoui tnicivshy 

Pla^Eie angle 
Sollii angle 

OUANTITV 

Force 

Gfkerey 
Powtr 

FItjA 

Fkyji doiaily 

Frequcii£^y 

ElcciriL; LDnductaiiOfl 
BlfiCtfomoLivc (oiM 



Unit 


Stmwl 




mcLft 


m 




klloBraiQ 


ka 




KCOOd 


t 




AEnper« 


A 




kelvjti 


K 




caDd«Ja 


cd 




otDle 


■mil 




Unit 


Symrol 




radian 


rad 




iteradi^n 


tr 




Unjt 


SVMBOt 


DfiFINlTldX 


djcwldd 


N 


1 xS - ] ks.tij/i* 


juUlfl 


J 


1 J ^1 N.ni 


watt 


W^ 


I W - 1 J^t 


wcbcjF 


Wb 


1 Wb - I V.s 


tcsla 


T 


IT =1 WfafftJ* 


herlft 


Hz 


1 Hz - 1 ci'»(4'') 


aicnieaa 


S 


IS - J A/V 


vol! 


V 


IV = 1 W^A 


paMial 


Pa 


I Pa ^ 1 N.'iD- 
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